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It is shown that general relativity coupled to nonlinear electrodynamics (NED) with the Lagrangian L(F) , F — 
F^vF^ having a correct weak field limit, leads to nontrivial spherically symmetric solutions with a globally regular 
metric if and only if the electric charge is zero and L(F) tends to a finite limit as F — > oo . Properties and examples 
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solutions are compared with their electric counterparts. A duality between solutions of different theories specified 
in two alternative formulations of NED (called the FP duality) is used as a tool for this comparison. 
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1. Introduction 

General relativity, despite its nonlinearity, is apparently 
lacking an effective self-restriction mechanism, and the 
existence of singularities seems to be its inevitable, though 
undesired, feature. Reasonable, regular solutions for 
macroscopic bodies such as stars are obtained with mat- 
ter whose pressure opposes gravity, whereas microob- 
jects, extreme states of matter and/or strong gravita- 
tional fields probably need a purely field description. 

The choice of a field source able to do the job is a sep- 
arate task, and, in particular, for spherically symmetric 
configurations there are quite a number of nonexistence 
theorems jl| . Non- Abelian gauge fields yield regular black 
hole solutions Q, but they are known only numerically. 
The regular black hole solution of Ref. || with a de Sitter 
core is expressed in terms of pressure and density rather 
than fields. An especially attractive class of field the- 
ories for seeking regular models is nonlinear electrody- 
namics (NED) with gauge-invariant Lagrangians L(F) , 
F = F iiV F lLV , since its energy-momentum tensor (EMT) 
T£ has the symmetry Tq = T± and is thus insensitive to 
boosts in the radial direction, which is a property of vac- 
uum j^, Q. Such theories, in particular, the Born-Infeld 
NED, recently gained much attention as limiting cases 
of certain models of string theory (see || for reviews). 
It has been shown, however |^], that in NED with 
any L(F) such that I ~ F at small F (the Maxwell 
weak- field limit), static, spherically symmetric general- 
relativistic configurations with a nonzero electric charge 
cannot have a regular center. As will be shown below, the 
same is true for dyonic configurations, where both electric 
and magnetic charges are present; regular solutions with 
wormholc topology also cannot exist for this system. 

The prohibition does not concern purely magnetic so- 
lutions, and, quite surprisingly, there is a whole class of 
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regular solutions with a nonzero magnetic charge. Regu- 
larity is here understood in a physical sense: although F 
is infinite at the center, the EMT is finite and the metric 
is regular (at least C 2 ) while the forces applied to test 
particles are finite everywhere and vanish at the center. 
The main aim of this paper is to present and to discuss 
these solutions. We will also compare them with their 
electric analogs, in particular, with the solutions found 
by Ayon-Beato and Garcia [8-10]. For this comparison 
we use a duality between spherically symmetric solutions 
of different NED specified in two alternative (F and P) 
frameworks: the original, Lagrangian framework and an- 
other one, obtained from it by a Legendre transformation 
& 

Throughout the paper all relevant functions are as- 
sumed to be sufficiently smooth, unless otherwise indi- 
cated. 

2. Nonexistence theorems 

Let us begin with a proof of two simple nonexistence the- 
orems extending the results of Q [7| . 

Consider NED in general relativity, with the action 

5 = lib / d ^ x ^~ 9 [R ~ L{F)l F = FfivF ^' (1) 

where R is the scalar curvature, F^ — d^A^ — d v A^ is 
the electromagnetic field, and L is an arbitrary function 
leading to the Maxwell theory at small F: L(F) w F as 
F —>■ 0. The tensor F^ v obeys the dynamic equations 
and the Bianchi identities, 

V^LpF^) = 0, V/F"" = 0, (2) 

where an asterisk denotes the Hodge dual and Lp — 
dL/dF. 
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A static, spherically symmetric metric can be written 
in the general form 

dg 2 = e 2j(u) dt 2 _ e 2 a (u) du 2 _ r 2^ dQ 2 (3) 

where dfl 2 = dO 2 + sin 2 8 dip 2 and u is a radial coordi- 
nate. A regular center, by definition, takes place where 
r = , if all algebraic curvature invariants are finite there 
and, in addition, there is a correct limiting value of the 
circumference to radius ratio (local flatness at the center, 
the absence of a conical singularity). If one chooses the 
curvature coordinates, u = r , the local flatness condition 
reads e°( 0) = 1 . 

The tensor F^ v compatible with spherical symmetry 
can involve only a radial electric field Fqi = — Fiq and a 
radial magnetic field F23 — — F32 ■ Eqs. give 

r 2 e a+ ^L F F 01 = q c , F 23 = q m sin (4) 

where q c and q m are the electric and magnetic charges, 
respectively. As follows from (||), 

f e d ^2F 01 F w = 2q 2 L- 2 r- i >0, (5) 

/ m d = f 2F 23 F 23 = 2g 2 1 r- 4 >0, (6) 

and the Einstein equations may be written in the form 

- Gl = T; = —2LpF llo ,F va + \%L (7) 
= i diag(L+2/ e i F , L+2f c L f , L-2f m L F , L-2f m L F ). 

(8) 

Theorem 1. The field system ^), with L(F) having a 
Maxwell asymptotic (L — > 0, Lp — > 1 as F — > 0), does 
not admit a static, spherically symmetric solution with a 
regular center and a nonzero electric charge. 

Proof. Since the Ricci tensor for the metric (^) is diag- 
onal, the invariant R^R* 1 " = R^Rt^ is a sum of squares, 
hence each component i?^ (no summing) is finite at a 
regular space-time point. Then each component of the 
EMT T£ is finite as well, hence, as follows from (0), 

(f e + f m )\L F \ <oo. (9) 

Suppose first that q m = 0, q c 7^ and thus f m = 
and F = — f e . Therefore, by (||) and (||), at a regular 
center FL F is finite whereas FL 2 F —> 00. Combined, 
these conditions lead to F — * and L F — * 00, that is, 
a strongly non-Maxwell behavior at small F , Thus for 
purely electric fields the theorem is valid. 

Suppose now q c ^ and q m =^ 0. Then (^|) should 
hold for f c and / m taken separately. As stated previ- 
ously, this condition applied to / c combined with (|5|) 
leads to L F — > 00 . But / m also tends to infinity as 
r — > 0, so even stronger f nl L F — * 00, violating (|^). The 
theorem is proved. 

A regular center is, however, not a necessary feature of 
a regular static, spherically symmetric space-time: there 



might be no center at all. Let us find out whether or not 
our system ([[J) can behave like this. 

If a center is lacking, the metric function r = {—g22) 1 ^ 2 
is restricted below by some minimum value r* > . There 
are then two opportunities: 

(i) r(u) has one or several minima (as it happens, e.g., in 

a wormhole) ; r itself then fails to be an admissible 
coordinate since it takes equal values at different 
spheres; 

(ii) r(u) monotonically approaches r m i n as u tends to a 

certain limiting value u* ; then the space-time can 
be nonsingular if it ends with a horn. 

By definition, a horn in a static, spherically symmet- 
ric space-time with the metric (^j is a region where, as u 
tends to some value u* , r(u) ^ const and gu = e 21 ^ 
have finite limits while the integral I = J e a du diverges. 
In other words, a horn is an infinitely long (3-dimensional) 
"tube" of finite radius, with the clock rate remaining fi- 
nite everywhere. It has an infinite spatial volume, and 
geodesies are infinitely continued along it as if in a worm- 
hole throat of unlimited length. The condition r ^ const 
discards purely cylindrical space-times, sometimes called 
flux tubes, which have no asymptotics. The definition 
of a horn follows the papers by Banks et al. [Q, where 
"horned particles" were discussed as possible remnants of 
black hole evaporation. 

The following theorem shows that these opportuni- 
ties can never be realized if the EMT has the "vacuum" 
property Tq = , and, in particular, for our system (|l|). 

Theorem 2. Let the metric ^) obey the Einstein equa- 
tions with an EMT satisfying the condition Tq = T^ . 
Then (i) the function r(u) cannot have a regular mini- 
mum and (ii) the space-time cannot contain a horn. 

Proof. Let us choose the u coordinate by fixing the 
condition a = lnr — 7. The ([[) — (J) Einstein equation 
takes the form d 2 (In r) / du 2 — . Therefore 

r(u) = e au+b , a,b = const, (10) 

so that r is either a constant (if a = 0) or a strictly 
monotonic function. This proves item (i). 

Suppose now that there is a horn. Then, by the above 
definition, a / 0, and a finite limit of r in Eq. ( [To| ) as 
u — > u* means that u* is also finite. On the other hand, 
since J e a du diverges as u — > u* , it follows that a —> 
00 , which, by our coordinate condition, can only occur if 
7 — * —00 , contrary to what was assumed. This completes 
the proof. 

Some remarks are now in order. First, the absence of 
wormhole solutions can also be proved from the known 
fact that a static wormhole throat always implies a vio- 
lation of the null energy condition [[[3) . This condition in 
our case reads Tq — T± > and is (marginally) observed 
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by the identical zero at the left-hand side. Our proof is, 
however, more direct and explicit. 

Second, the opportunity 7 — > —00, mentioned in the 
proof of item (ii) of Theorem 2, generically corresponds to 
an event horizon. [It can happen in principle that e 7 — > 
at finite r does not imply a horizon. Our system ([!]) with 
smooth L{F) does not admit such cases, while horizons 
do exist, as can be seen from the solution below.] 

Third, the above theorems do not contain an asymp- 
totic flatness requirement, the proofs being of local na- 
ture. Therefore both theorems are readily extended to 
general relativity with a cosmological constant (the con- 
dition Tq = Tl holds in its presence as well), where the 
spatial asymptotic can be de Sitter or anti-de Sitter. 



3. Regular magnetic solutions 

Both theorems were proved without entirely solving the 
Einstein equations. For our system (|l|), however, an ex- 
act solution can be obtained by quadratures in the general 
spherically symmetric case [||. Indeed, the Maxwell-like 
equations are already integrated. Let us choose the curva- 
ture coordinates u = r , which is now safe since we know 
that r has no extrema. Due to Tq = Tl , the correspond- 
ing Einstein equation gives d(a+^)/dr = a+7 = 
for a proper choice of the time unit. It remains to write 
the well-known relation for a(r) in terms of the energy 
density Tq , which follows from the ([J) Einstein equation: 



=27 
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Air) 



2M(r) 



M(r) = 



- /T °(r)r 2 dr. 



(11) 



Possible horizons occur at zeros of A(r) . 

The only nontrivial case not covered by the theorems 
is q c = , g m 7^ , when there is still a hope to obtain a 
regular center. In this case the metric has the form (||) 
with ( |ll| ) where 



M(r) 



\ J L(F)r 2 dr 



(12) 



and F = 2q 2 ri /r i . It is easily seen that a space-time with 
a regular center is indeed obtained for any L(F) such that 
L — > Too < 00 as F — > 00 when one integrates in Eq. jl^ ) 
from to r . Integration from to 00 then gives a unique 
mass m = M(oo) providing a regular center for given 
q m ; hence the entire mass is of electromagnetic origin. A 
finite value of M(oo) is guaranteed by our assumption of 
a Maxwell asymptotic of L(F) at small F: as r — > 00, 
L ps 2g 2 1 /r 4 , and the integral ( |l2| ) converges. 

The EMT in (@) near r = takes the form T£ = 
5<5^Too(l + o(l)), and the metric is approximately de 
Sitter (A(r) = 1 — Ar 2 /3 + o(r 2 ) ), with the cosmological 
constant A = Loo/2. The Riemann tensor at r — coin- 
cides with that of de Sitter space, so one need not explic- 
itly calculate the curvature invariants to prove that the 



space-time is regular. The metric is at least C 2 smooth 
at r = but, depending on L(F) , may be even C°° 
smooth, as will be seen in the example (p9[). 

Suppose that L(F) and the mass have been chosen 
as described. The space-time is then globally regular 
and can include horizons corresponding to zeros of A(r) , 
whose number and character determine the global struc- 
ture [note that A(0) — A(oo) = 1]. In the absence of 
zeros, there is a regular Dirac-type magnetic monopole 
solution. The occurence of two simple zeros leads to 
the conventional Reissner-Nordstrom black hole struc- 
ture, with the singularity replaced by a regular center. 
The intermediate case of one double zero gives the ex- 
tremal Reissner-Nordstrom structure. Models with more 
numerous horizons can be constructed as well. 

Due to the above theorems, our magnetic black holes 
and monopoles are the only types of static, spherically 
symmetric solutions to with a globally regular metric. 

The only trouble is an infinite magnetic induction B 
(B 2 = / m /2 = F/2) at r = 0, whereas the magnetic field 
intensity H , obtained as a generalized momentum from 
the Lagrangian, is well-behaved everywhere including the 
center: H 2 = f m L 2 F /2 — ► as r — * 0. To judge whether 
or not the center is regular from a physical viewpoint, one 
should estimate the force applied to a charged test par- 
ticle moving in the field under consideration. This test 
charge may be electric or magnetic since both are admit- 
ted by our assumptions. In a consistent approach, the 
equations of motion for a test particle (as well as those 
for an extended body) in nonlinear field theory should fol- 
low from the field equations and may be deduced along 
the lines of Refs. |L4], ^aj. Namely, the 4-force vector is 
found as the total momentum flow J T^n^dS through a 
closed surface surrounding the particle, where n M is the 
unit normal to such a surface and is the total EMT 
of the summed electromagnetic field of the background 
static configuration and the test particle itself. An es- 
timate in a proper approximation, taking into account 
the weakness of the particle's field, shows that this force 
is finite everywhere and vanishes at r = . Therefore, 
despite the divergent B , our magnetic solutions may be 
called globally regular. 

4. FP duality and electric solutions 

Let us consider for comparison the electric analogs of our 
magnetic solutions. This is of particular interest since 
recently Ayon-Beato and Garcia [8-10] suggested some 
examples of such solutions, describing configurations with 
q c 7^ 0, q ln = and a regular center. The properties of 
these solutions evidently contradict Theorem 1, but they 
only seem to circumvent it since, as we shall see, any 
model like those of [8-10] needs different Lagrangians in 
different ranges of the radial coordinate and therefore fails 
to be a solution for a particular Lagrangian L(F) . 

The solutions of [8-10] were found using an alter- 
native form of NED (to be called the P framework), 
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obtained from the original one (the F framework) by 
a Legendre transformation: one introduces the tensor 
Pfiv = LpF^ with its invariant P = P iLV P^ v and con- 
siders the Hamiltonian-like quantity H = 2FL p — L as a 
function of P; the theory is then reformulated in terms 
of P and is specified by H(P) 11 . One then has: 



L = 2PH P -H, L F H P = 1, F = PH\ 



(13) 



with Up = dH/dP. Eqs. (§) and the EMT (@) are rewrit- 
ten in the form 



V lL P^ = 0, 



V M (W P *FH = 0, 



(14) 



T" = ~\ diag (H-2p m H P , n-2 Pm n P , 



H+2 Pc Hp,H+2 Pc Hp) (15) 
where, by (|l4|), in the spherically symmetric case, 

p c - 2P m P w = 2q 2 Jr\ (16) 



p m = 2P 23 P ZA = 2qiJip/r 



(17) 



Comparing Eqs. (§, (|)-(|) with @-(Jl7|), one sees 
that they coincide up to the substitutions 



{5p 



F F 



L(F)} 



-H(P)}. 
(18) 

In other words, there is a duality between spherically sym- 
metric solutions written in the F and P frameworks: any 
solution for a given Lagrangian L(F) , characterized by 
a certain metric function A(r) and the field components 
Fqi and F23 , has a counterpart with the same A{r) but 
with F substituted by — P , L by — TC, F i by P23 and 
F23 by Pqi , and conversely. The functional dependence 
—Ti.(—P) in the dual solution is the same as L(F) in 
the original solution, but the choice of the function L(F) 
itself is not restricted. 

(It should be stressed that this duality, to be called 
FP duality, connects solutions of different theories: 
given L(F), the functional dependence H-(P) = 2FLp — 
L is in general quite different from L(F) , an evident 
exception being the Maxwell theory, where L = F = 
H = P and the present duality turns into the conven- 
tional electric-magnetic duality. So the FP duality has 
nothing to do with the electric-magnetic one studied in 
16 1, where the field equations of a specific theory 



Refs. 

were required to be duality invariant, and this condition 
selected a narrow class of Lagrangians.) 

In particular, any regular magnetic solution obtained 
for given L(F) has a purely electric counterpart with a 
similar (up to the sign) dependence H(P) . The metric 
has the form ( |IT|) with 



M(r) 



\ J H{P)r 2 dr. 



(19) 



Given H{P) , one should substitute P = -2q 2 /r 4 . A 
regular center exists if and only if H. has a finite limit as 
P — > — 00, and a mass that provides regularity for given 




Figure 1: An example of qualitative behavior of H(P) , F(P) , 
and L(F) in an electric solution 



q c is found by integration in (|l^) from to 00 . This is 
how the regular solutions of [8-10] were obtained with the 
following choices of Tt(P) : 



H(P) 



P 



1 - 311 



6 



n 



5/2 



(20) 



(1 + n) 3 9 2 sVi + n 

§:W(P) = F/cosh 2 (sVn), (21) 



where II = ^/—q 2 P/2 and s = \q\/(2m), q — q c and 
m being free parameters identified with the charge and 
mass of the configuration, respectively. The functions 
(^ti|)-(^2|) behave like P at small P , tend to finite limits 
as P — > —00 and thus lead to regular metrics. 

Yet the P framework is secondary: the Lagrangian 
dynamics is specified in the F framework, and, since the 
F i— > P transition is a mere substitution in the field equa- 
tions, the two frameworks are only equivalent where the 
function P(F) ismonotonic. Recalling the proof of Theo- 
rem 1 for q c 7^ 0, q m = 0, one sees, however, that for any 
regular solution with a Reissner-Nordstrom asymptotic 
the function F(P) — — / < vanishes at both r = 
and r — > 00 . So it inevitably has at least one minimum 
at some P = P* < . It can be shown directly that at an 
extremum of F(P) where F = F* < the derivative Lp 
has the same finite limit as P — > P* + and P — > P* — , 
while Lpp tends to infinities of opposite signs. Therefore 
the function L{F) suffers branching, and its graph forms 
a cusp at F — F* ; different functions L(F) correspond 
to P > P* and P <P*. 

Another kind of branching occurs at extrema of Ti.(P) , 
if any. There F(P) behaves generically as (P — P*) 2 
while Lp — > 00, and a graph of L{F) smoothly touches 
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the vertical axis F = . The number of Lagrangians on 
the way from infinity to the center equals the number of 
monotonicity ranges of F(P) . 

All this is readily seen for specific examples. A quali- 
tative picture for the choice ( |2l| ) is shown in Fig. 1. 

In the simplest case when H(P) is monotonic (e.g., 
like tanhP), L(F) has only two branches OP\ and 
P1P2 , and P2 already corresponds to r = 0. 

Thus any regular electric solution, being well-behaved 
with respect to the field equations in the P framework, 
corresponds to different Lagrangians in different parts of 
space. This problem is absent for magnetic solutions since 
they are obtained directly in terms of L(F) . 

5. Comparison of effective metrics 

The troubles with the electric solutions concern only the 
properties of NED, while the metric is quite well-behaved. 
The same is true for the electric field Fo\ . However, ter- 
mination of a theory with given L(F) implies some vio- 
lent electromagnetic phenomena. For their understanding 
let us consider the effective metric introduced by Novello 
et al. 0, H| 



ft"" = .9cff = 9^L F - AL FF F» a F° 



(23) 



As shown in |L7| [l8| , NED photons propagate along null 
geodesies of this metric. For the space-time metric (11), 
with a purely electric field the effective metric reads 



ds. 



huvdx^dx" = — 



A(r)dt 2 - 



dr 2 



!ff - $ 

$ = L f + 2FL ff =Hp/F p . 



Mr) 



Lp 



(24) 



At an extremum P = P* of F{P) where F^O (in par- 
ticular, at the inevitable first minimum) one has $ — ► 
since Fp — > while Hp is finite. This leads to a curva- 
ture singularity of the effective metric, at least if P* is 
not located on a horizon, A ^ . Another kind of singu- 
larity of the metric (|24|) accompanies possible extrema of 
Ti.(P) . All this is verified by calculating the Kretschmann 
scalar K . Even more importantly, according to |l8| , if 
a NED photon comes from an emitter at rest at point 
1 to an observer at rest at point 2, the corresponding 
frequencies f\ and fa are related by 



h 



Ml) / M2) _ *(2)/v^(2) 



v^(T)/ v^(2) mivw) 



(25) 



where the second equality corresponds to the metric 
(|2^). If $(2) = 00 [as it happens at a termination 
point of L(F)], then photons coming there are infinitely 
blucshifted and one may expect that they eventually lead 
to a real space-time singularity. 

For a magnetic solution, instead of (p3), we get 



A{r)dt 2 



dr 2 

W) 



—dVL 2 (26) 



where again $ = Lp + 2FL F F . Instead of fl2 



k 
ft 



Lp{2)/^MT) 



L F {l)/^W) 



(27) 



At the center (r = 0, A = 1) both L F and $ vanish, the 
coefficient /122 — > 00, i.e., behaves as if in a wormholc, 
whereas /loo — ► 00 , which means that photons arriving 
there, if any, would be infinitely redshifted [see Eq. (p7|)]. 
Actually photons cannot reach a place where Lp = , as 
can be seen from an integral of their geodesic equation 



L 



-2-2 



2l;2 



[A(r)$y]l 



(28) 



where the overdot is a derivative in the affine parameter, 
and e and / are the photon's constants of motion char- 
acterizing its initial energy and angular momentum. All 
curvature invariants of the metric ( p6| ) vanish at r = . 
It is indeed a perfectly quiet place despite an infinite F . 

Some peculiarities, however, occur on the way from 
infinity to the center. There is always a sphere r — r* on 
which $ = . It can be seen as follows: $ may be rep- 
resented as $ = 2y/F(\/FLp) F . The quantity \/FL F 
vanishes at both r = and r = 00 and is nonzero be- 
tween them, hence it has at least one extremum at F =/= 
— this is where $ = . The metric (^6|) , due to blowing- 
up of the coordinate spheres, has a singularity there, but 
the latter is actually unnoticed by NED photons, as is 
evident from Eq. (Es) . Generically L F ^ where $ = , 
therefore the photon frequency also remains finite. The 
meaning of the very fact of a curvature singularity of the 
effective metric is yet to be understood. 

If L F = at some F > 0, this also causes a singu- 
larity of Eq. ( |26| ) , which acts as a potential wall (mirror) 
for NED photons, as is seen from (p8|). Accordingly, ( 570 
shows that they are infinitely redshifted: fa vanishes if 
L F (2) = 0. No photons from outside can thus approach 
the center. 

All this is in striking contrast to the picture obtained 
for an electric source — we now have potential walls in- 
stead of wells and redshifts instead of blueshifts. 

6. Example 

To have a specific example of a regular magnetic solu- 
tion, let us employ the above FP duality and consider, 
with slight modifications, the dependence (pl|), substi- 
tuting — H with L and — P with F . An advantage of 
Eq. @ [as well as @ and @] is that it leads to a 
closed form of M(r) and A(r) . Let us, however, slightly 
modify it, excluding an explicit dependence of L on m 
and q: they should be integration constants, while L may 
only contain fundamental constants or those originating 
from a deeper underlying theory. Moreover, to be able to 
describe systems with both electric and magnetic fields, 
where F (and P) can have both signs, let us replace — P 
with \F\ rather than F . So we put 



L(F) = F/cosh 2 (a|F/2| 



l/4\ 



const. 



(29) 
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The use of |F| violates analyticity of L at F — 0: as 
required, Lp(0) = 1, but Lpp contains the discontinu- 
ous term — a 2 |2F| -1 / 2 signF. Though, in the range of 
interest, F > 0, this L(F) is well-behaved. Integration 
in Eq. (|l^) gives for a regular solution 



M(r) 



13/2 



2ft 



1 — tanh ■ 



(30) 



so that to = M(oo) = \q\ 6 ' /2a (q = <7 m )> and some 
relations from are formally restored. In particular, 
the minimum value of A(r) = 1 — 2M(r)/r [recall that 
A(0) = 1] depends on the ratio £ = m/\q\, so that A m ; n 
is negative for £ > £ ~ 0-96 (we deal with a black hole 
with two horizons), zero for £ = £o (an extremal black 
hole with one double horizon) and positive for £ < £ 
(a regular particle- like system). It is of interest that, 
given any specific value of the constant a in Eq. ( f29| ) , we 
can obtain all three types of solutions depending on the 
charge value: we have a non-extremal or extremal black 
hole if \q\ < 4a 2 /£q, or we have a particlelike solution 
(a monopole) otherwise. Despite the restriction imposed 
by the regularity condition, one finds all three types of 
regular solutions. This feature seems to be quite generic 
for proper nonlinear Lagrangians. One can also verify 
that the properties of the effective metric ( |26| ) confirm 
the above general observations. 

One can also notice that, due to an exponential decay 
of M(r) in Eq. (|30|), the metric is in this case C°° smooth 
at r = . 

7. Concluding remarks 

A more complete description of the properties of the 
present regular NED solutions, as well as others, re- 
quires a better understanding of the long-standing and 
non-trivial problem of motion of charged bodies in NED, 
probably following the lines of Refs. (lj, [ll| [is|] . 

One more subject of interest for further study is 
the inclusion of another electromagnetic field invariant, 
*F tJ - u F fll/ , into the Lagrangian in addition to F. This 
invariant is involved, in particular, in the Born-Infeld 
and Heisenberg-Euler NED Lagrangians; its appearance 
should be able to widen the diversity of regular black 
hole and monopole solutions. Related subjects are the 
FP duality between solutions of different theories in- 
volving both invariants and a possible extension of this 
duality to non-spherically symmetric configurations. 
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